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ABSTRACT 

We report on the first detection of 13 C enhancement in two B[e] supergiants in the Large 
Magellanic Cloud. Stellar evolution models predict the surface abundance in 13 C to strongly 
increase during main-sequence and post-main sequence evolution of massive stars. However, 
direct identification of chemically processed material on the surface of B[e] supergiants is 
hampered by their dense, disk-forming winds, hiding the stars. Recent theoretical computa- 
tions predict the detectability of enhanced 13 C via the molecular emission in 13 CO arising 
in the circumstellar disks of B[e] supergiants. To test this potential method and to unam- 
biguously identify a post-main sequence B[e]SG by its 13 CO emission, we have obtained 
high-quality A"-band spectra of two known B[e] supergiants in the Large Magellanic Cloud, 
using the Very Large Telescope's Spectrograph for INtegral Field Observation in the Near- 
Infrared (VLT/SINFONI). Both stars clearly show the 13 CO band emission, whose strength 
implies a strong enhancement of 13 C, in agreement with theoretical predictions. This first ever 
direct confirmation of the evolved nature of B[e] supergiants thus paves the way to the first 
identification of a Galactic B[e] supergiant. 
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1 INTRODUCTION 

B[e] supergiants (B[e]SGs) are massive and luminous B stars 
with strong, non-spherical, disk-forming winds. The conditions 
in terms of density and temperature within these disks are ideal 
for efficient molecule and dust formation, traced by CO band 
emission (McGregor et al. 1988; Morris et al. 1996) and a strong 
near and mid-infrared excess Cickgraf et al.ll 1986T) . However, the 
usual distance issue renders it difficult to obtain useful luminos- 
ity es timates to classify a Galactic B[e ] star definitely as B[e]SG 
(e.g.. lLamers et alj[l998l ; Ikraus 2009). An additional complica- 
tion arises from the position of Galactic B[e]SG candidates in 
the empirical Hertzsprung-Russell diagram, which they share with 
the pre-main sequence Herbig Ae/Be stars. Even worse, the lat- 
ter are also surrounded by gas and dust disks JWatersI 120061 : 
IWaters & W aelkens 1998), giving rise to identical observable fea- 
tures such as the infrared excess and the well-pronounced C O band 
emission JBik & Thill2004lThi et al]2005l : lBik etai]|2006tl . For an 
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unambiguous classification of a Galactic B[e] star as B[e]SG, a re- 
liable tracer based on chemical processing during the evolution of 
a massive star is thus needed. 

Based on stel l ar evolution models f o r Gal a ctic stars 
dSchaller et ai]|l992l ; iMevnet & Maederl I2005T) . iKrauj J2009h re- 
cently pointed out the significant enrichment of the stellar surface 
abundance in 13 C already during the main-sequence evolution of 
massive stars. However, while the surface of B[e]SGs is usually 
hidden within their dense winds and disks, IKrausI ( 120091) proposed 
that this enrichment should be detectable in the ejected material 
via strongly enhanced 13 CO band emission from the disk-forming 
winds of B[e]SGs. If true, we would have, for the first time, an easy 
and robust method to unambiguously distinguish between pre- and 
post-main sequence evolution and to classify Galactic B[e] stars 
as supergiants, because Herbig Ae/Be stars do not show any C 
enhancement. However, before this method can be applied to the 
Galactic B[e] stars, its validity needs to be tested on well-known 
B[e]SGs. The best-known sample of confirmed B[e]SGs is lo cated 
in the Magellanic Clouds JLamers et al.lll998l : lzickgraj|2006l) . 

To test the 13 CO method, we have observed two B[e]SGs in 
the Large Magellanic Cloud (LMC) which are known to display 
12 CO band emission originating in their circumstellar disks. In this 
Letter we report on the first definite detection of 13 C enrichment 
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Table 1. Parameters of the observed B[e]SGs. 



Object 
LHA 


Sp. Type 


T c s 
[10 3 K] 


[10 5 L Q ] 


[RqI 


M ini 
[M ] 


£(B - V) 


[kins- 1 ] 


inclination 


Reference 


120-S 12 
120-S 73 


B1I 
B8I 


23 
12 


2.2 
3.0 


30 

125 


~25 
25-30 


0.25 
0.12 


285 ± 10 
264 ±3 


intermediate 
~ pole-on 


Zickaraf et al. (1986) 
Stahl et al. ( 19831 
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Figure 1. Parts of the observed, flux-calibrated SINFONI iC-band spectra 
of the B[e]SG stars S 73 (top) and S 12 (bottom) covering the Br7 line and 
the first-overtone bands in both 12 CO and 13 CO. The wavelengths of the 
band heads are indicated by vertical ticks. 



from B[e]SGs via enhanced CO band emission whose presence 
is direct evidence for the evolved nature of our programme stars. 



2 OBSERVATION AND REDUCTION 

Between October and November 2009, we have obtained high- 
quality if-band spectra (1.95 - 2.45 /im) of the two LMC B[e]SG 
stars LHA 120-S 12 and LHA 120-S 73 (in the following S 12 and 
S 73 respectively), using the Spectrograph for INtegral Field Ob- 
servation in the Ne ar-Infrared (SINFONI) (Eisenhaue r et alj|2003l : 
iBonnet et alJuOOot) . Observations were carried out in an AB pat- 
tern with a field of view of 8" x 8". B-type standard stars for flux 
calibration were observed at similar airmass. 

Data reduction was performed with the SINFONI pipeline 
(version 2.0.5). Raw frames were corrected for bad pixels, distor- 
tion, and flatfield, followed by the wavelength calibration. For flux 
calibration, we extracted the stand ard-star spectra and scaled them 
with the according K urucz mode l (Kurucz 1993) to their 2MASS 
i\~ s -band magnitude dSkrutskie et al.l 12006) to obtain calibration 
curves. The final flux-calibrated spectra have a signal-to-noise ratio 
S/N ~ 500 and a spectral resolution R ~ 4500 at 2.3 /im. Fig- 
ureQ] shows parts of the if-band spectra covering the Br7 line of 
hydrogen and the CO bands. 
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Figure 2. Change in 12 C/ 13 C surface abundance ratio during the evolu- 
tion of rotating massive stars at LMC metallicity, based on the models of 
iMevnet & MaedeJ 12005]) . 



3 ANALYSIS AND RESULTS 

The two LMC B[e]SGs listed in TableQ]were selec ted because they 
are known to display strong 12 CO band emission l IMcGregor et all 
1988). The pole-on or intermediate orientation of their disks guar- 
antees that we see the complete CO emitting disk area. The spectral 
range of SINFONI covers the first four band heads of 12 CO and the 
first two of 13 CO. Visual inspection of the spectra (Fig.QJ reveals 
that the 1 CO band heads are prominent in both stars; while 1 CO 
is clearly visible in S 73 (top panel), its presence in S 12 is slightly 
less obvious. 

We corrected the spectra for the systemic velocity and dered- 
dened them with the E(B — V) values from TableQJapplying the in- 
terstellar extinction curve of Howarth (1983) and using Av/E(B — 
V) = 3.1. To obtain a flux-calibrated, quasi-pure CO emission 
spectrum, we fitted and subtracted a linear pseudo-continuum long- 
wards of 2.25 fim. The spectra of the two stars do not show indica- 
tion for significant emission of the hydrogen Pfund series, but some 
minor contribution cannot be excluded and might have to be taken 
into account later on in the analysis. 

Inspection of the CO spectra reveals that the band heads of 
12 CO all have about equal strength, which means that, surprisingly, 
the CO gas is rather cool (T C o < 3000 K). We will come to this 
aspect again in Sect.|4] On the other hand, the excitation of pro- 
nounced band head structures requires CO temperatures in excess 
of ~ 2000 K. Therefore, the location of the CO emitting region 
must be closer to the star than the dusty regions, because dust sub- 
limates at ~ 1600 K (e.g., Lodders 2003). Moreover, the rise of 
the CO quasi-continuum towards longer wavelengths indicates sub- 
stantial optical thickness of the emission. To model the CO bands 
we can thus assume that the emission is in local thermodynamical 
equilibrium (LTE) and, preempting results from Sect.|4] that the CO 
gas is located in a ring around the star which has constant tempera- 
ture and column density. These assumptions are reasonable since 
the hottest CO component, which arises close to the star where 
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Figure 3. Top.' Synthetic spectra of the CO bands (red) and of the Pfund series (blue). Bottom: Combined CO band plus Pfund spectra (red) overlayed on the 
observed (black), flux-calibrated CO band spectra of S 73 (left) and S 12 (right). The vertical ticks in the respective panels indicate the positions of the 13 CO 
band heads. 



usually the disk density is highest, dominates the total CO spec- 
trum completely. F or d etails on the C O band modeling we refer to 
iKraus et al.l feOOd) and lKrausI d2009h . 

To constrain the range of the 12 C/ 13 C abundance ratio in the 
disk, we took expected 12 C/ 1 3 C surface a bundance ratios from 
the stellar evolution models of iMevnet & Maedea 1 120051) for mas- 
sive stars with LMC metallicity and initial rotation velocities of 
300 km s -1 . Fig.|2] shows how this ratio decreases with time for 
stars of different initial mass. The stars S 12 and S 73 have evolved 
from progenitors with Mini = 25 ... 30 Mq (see TableflJ whose 
surface abundance should have dropped to a value 12 C/ 13 C < 20 at 
the turn-off from the main sequence and decreases down to a value 
of ~ 4 at late evolutionary phases. We computed synthetic CO band 
spectra for temperatures in the range 2000 - 3000 K, 12 C/ 13 C ra- 
tios varying from 4 to 20, and CO column densities higher than 
10 21 cm -2 to account for the optical depth effects. We did not find 
any indication for additional line broadening from either Keple- 
rian rotation or outflow signatures, meaning that the CO gas has no 
significant line-of-sight velocity component, in agreement with an 
(almost) pole-on orientation of the disks. Therefore, we used the 
thermal velocity as the intrinsic gas velocity and convolved the re- 
sulting spectrum with the resolution of SINFONI. Pure CO models 
that match the observations best are shown in red in the top panels 
of Fig.[3] with the parameters listed in Table[2] 

The presence of some additional emission features, especially 



at the red end of the spectrum, can be explained by the (if weak) 
recombinatio n lines of the hydr ogen Pfund series. Following the 
description of lKraus et alj 1 12000]) . we modeled the Pfund series as- 
suming a constant electron density and an electron temperature of 
T e = 10 4 K. The lack of Pfund lines higher than Pf 34 (S 73) and 
Pf 31 (S 12) means that due to pressure ionization in a high-density 
gas, recombination into such high-energy states is suppressed. This 
is in good agreement with the electron densities of n e > 10 9 cm -3 
needed for the lines to be formed in LTE conditions, as is obvious 
from their intensity ratios. As in the case of the CO bands, no addi- 
tional line broadening is required to match the width of the Pfund 
lines. This could indicate that the Pfund emission is generated ei- 
ther in the disk but much closer to the star than the CO bands, or 
in the polar wind regions so close to the star that the wind velocity 
is still very small compared to the velocity resolution of SINFONI. 
The resulting synthetic Pfund series are included as the blue lines 
in the top panels of Fig.[3] and the combined model spectra of CO 
bands and Pfund lines are overplotted in red on the observed spectra 
in black in the bottom panels. 



4 DISCUSSION 

Our analysis has shown that the spectra of the B[e]SGs S73 and 

5 12 display emission from the 13 CO band heads (see Fig.[3} which 



© 2010 RAS, MNRAS 000,[[H5] 



L4 A. Liermann et al. 



Table 2. Best-fit parameters for the CO band emission from the disks of 
S 12 and S 73. 



Object 
LHA 



Too -?Vco 

[K] [10 21 cm" 2 ] 



12 C/ 13 C Aco cos i 

[AU 2 ] 



120-S 12 
120-S 73 



2800 ± 500 
2800 ± 500 



2.5 ±0.5 
3.5 ±0.5 



20 ±2 
9±1 



2.58 ±0.15 
21.0 ±0.3 



means that their circumstellar disk material is strongly enriched 
with 13 C; hence, these stars are indeed evolved post-main sequence 
objects. 

Our model computations have also confirmed the relatively 
low temperature of the CO gas as indicated by the about equal 
strengths of the 12 CO band heads. Such low temperatures are not 
expected within the canonical picture of B[e]SGs, which assumes 
that their disks are formed by enhanced mass flux from the stars' 
equatorial regions. If the disks of the two B[e]SGs studied here 
were continuously supplied, then the CO in the gaseous parts of 
the disk should be visible from the hottest component having tem- 
peratures close to the CO dissociation temperature of 5000 K. In 
this case the observable CO band spectrum, which always reflects 
the hottest available CO component, would be dominated by the 
peaks of the higher vibrational transitions occurring at longer wave- 
lengths, while t he lowest ban d head (at 2.3 ^im) would be strongly 
suppressed (see iKraus 2009). The apparent deficiency in CO gas 
hotter that ~ 2800 K therefore suggests that there is no CO gas 
closer to the star. Hence, the disks seen around S 12 and S73 can- 
not extend down to the stellar surface. Instead it would be more 
appropriate to assume that the stars are surrounded by a dense 
and cool ring of material. Support for such a scenario comes from 
the recent detection of a detached ring in quasi-Keplerian rotation 
arou nd the Small Mage llanic Cloud (SMC) B[e]SG star LHA 115- 
S65 JKrausetai]|2010|) . 

Alternative ly, gravitational darkening according to the von 
Zeipel theorem Jvon Zeipelll 19240 in rapidly rotating stars can lead 
to equatorial temperatures being sufficiently low (i.e., < 3000 K) 
for the hottest CO gas to not exceed ~ 2800 K. In this picture, 
the disk still could extend from the stellar surface out to large dis- 
tances, with the CO emitting gas located closest to the star within 
the innermost, hottest disk region. While this scenario is somewhat 
speculative (and heating of the disk through light from hotter re- 
gions of the star would have to be considered as well), we can es- 
timate the minimum rotation speed of the stars required to achieve 
an equatorial surface temperature of 3000 K. 

Stellar rotation leads the poles to be hottest. Assuming that 
both stars are seen pole-on and that their effective temperatures, 
T e ff, listed in TableQ]correspond to the polar values, requires a de- 
crease in T c fi from polar to equatorial values by factors of 7.7 (S 12) 
and 4 (S 73). This, in turn, requires the stars to rotat e with 3> 99% 
of their critical velocity (see Fig. 5 of |KrausJ|200g). The assumed 
pole-on orientation only provides us with lower limits of the real 
rotation velocity, but even in the unexpected case that both stars are 
seen perfectly pole-on, rotation rates close to or even at the critical 
limit are highly unlikely. For comparison, the two most rapidly ro- 
tating B[e]SGs known, the SMC stars LHA 1 15-S 2 3 and LHA. 115- 
S 65, rotate with > 75% of their critical velocity (Zickgraf 2000; 
IKraus et al. 2008, 2010). If our program stars rotated at ~ 75% of 
their critical velocity, resulting equatorial temperatures would still 
be ~ 72% of the polar values (cf. Table[TJ, i.e. much hotter than 
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Figure 4. Outer radius of the CO emitting ring as a function of inner radius 
(solid lines) for S73 (top) and S 12 (bottom). The maximum outer radius 
is set in both cases by the distance for dust sublimation (dashed lines). For 
S 12 two inclination angles for the disk are considered. 



3000 K. The deficiency of CO gas hotter than 2800 K is thus real, 
and the most plausible scenario is hence that the molecular gas is 
located within a detached ring of material rather than in a disk ex- 
tending down to the stellar surface. 

In Table[2] we include the sizes of the CO emitting regions 
(projected to the line of sight, Aqo cos i) as obtained from our best- 
matching models. From these sizes we estimate the radial extension 
of the emitting CO rings. We compute the outer radius of the ring as 
a function of its inner radius. The results for both stars are depicted 
by the solid lines in Fig.|4]indicating that for increasing inner radii 
the emission area converges towards infinitesimally narrow rings. 
With the stellar luminosities given in TableQ]and a n assumed dus t 
sublimation temperature of T su b — 1600 K (e.g., Loddersl l2003l) 
we find individual sublimation radii of i? S ub — 24 R* and ~ 86 R„ 
for S 73 and S 12, respectively. These maximum outer edges of the 
CO rings are shown in Fig.|4]as dashed lines. An estimate for the 
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inner edge of the CO ring is less precise. The CO gas has to have 
a temperature of about 2800 K, only slightly higher than the dust 
sublimation temperature but definitely much lower than the effec- 
tive temperatures of the stars. A reasonable estimate could thus be 
-Riii •> 10 Rt. for S 73 and (much) larger for S 12. For the latter we 
plotted the relation between inner and outer radius for two different 
values of disk inclination. However, independent of the exact value 
of Ri n (and i as in the case of S 12) we can conclude that the CO 
emitting area is concentrated within a rather narrow ring at short 
distance from the central star. 



5 CONCLUSIONS 

Using high-quality Tf-band spectra obtained with SINFONI, we 
detected, for the first time, a 13C enhancement in two B[e]SGs 
in the LMC, thereby confirming their evolved nature. The enrich- 
ment of their disk-forming wind material with 13 C and the result ing 
strong 13 CO band emission, as was suggested by Kraus (2 0091) . is 
thus shown to be a robust method for the identification of evolved 
stars. This method can now be applied to the Galactic B[e]SG can- 
didates to disentangle the real Galactic B[e]SGs from the more 
abundant Herbig Ae/Be stars. 

In our analysis, we also found a deficiency of hot CO gas in 
the disks of both stars. This can be explained in two ways. Ei- 
ther both stars rotate at their critical limit so that equatorial sur- 
face temperatures (and hence, maximum disk temperatures) drop 
below ~ 3000 K and no hot CO gas is expected, or the CO emit- 
ting region is located in a high-density, detached ring of material 
rather than in an equatorially enhanced, disk-forming wind. While 
the former scenario cannot be excluded, we consider the latter as 
the more realistic one since it also agrees with the recently found 
detached gas ring around the SMC B[e]SG star LHA 1 15-S 65 by 
iKrausetalJ 12010) . 
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